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O ■ ABSTRACT 
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We present a measurement of the D/H abundance ratio in a metal-poor 
^ \ damped Lya (DLA) system along the sightline of QSO SDSS1558-0031. The 

DLA system is at redshift z = 2.70262, has a neutral column density of log Nm = 
20.67±0.05 cm -2 , and a gas-phase metallicity [0/H]= —1.49 which indicates that 
deuterium astration is negligible. Deuterium absorption is observed in multiple 
Lyman series with a column density of log Ad = 16.19 ± 0.04 cm -2 , best con- 
q | strained by the deuterium Lyman-11 line. We measure log D/H = —4.48 ± 0.06, 

which when combined with previous measurements along QSO sightlines gives a 
O ' best estimate of log D/H = —4.55 ±0.04 where the la error estimate comes from 

O 



a jackknife analysis of the weighted means. Using the framework of standard 
big bang nucleosynthesis, this value of D/H translates into a baryon density of 
Qbh 2 = 0.0213 ± 0.0013 ± 0.0004, where the error terms represent the la errors 
from D/H and the uncertainties in the nuclear reaction rates respectively. Com- 
bining our new measurement with previous measurements of D/H, we no longer 
find compelling evidence for a trend of D/H with Ahi- 

> '■ 
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1. INTRODUCTION 

For the last decade, measurements of the primordial D/H ratio in QSO sightlines have 
provided increasingly more precise constraints on the cosmological baryon density. Although 
the measurement of D/H is simple in principle, compared to the other light elements pro- 
duced during big bang nucleosynthesis, finding those QSO absorption lines systems which 
are suitable for measuring D/H has proven observationally challenging (Tytler et al. 2000). 

For a QSO absorption system to show D/H, a number of criteria must be met (see 
Kirkman et al. (2003), hereafter K03, for a more detailed discussion). First, the hydrogen 
column density must be large enough (since D/H is of the order of one part in 10 5 ) such that 
deuterium can be observed using modern high-resolution spectrographs. Second, the velocity 
structure of the hydrogen absorption must be simple enough, ideally a single component of 
gas, so that the deuterium absorption is well resolved given the small 82 km s _1 offset from 
the hydrogen Lyman lines. Third, there can be little to no interloping Lja forest or metal 
lines at the position of the deuterium absorption, since such absorption strongly complicates 
attempts to constrain the deuterium column density. Unfortunately, Lja forest absorption 
is both ubiquitous and stochastic in high redshift QSO spectra. Finally, the background 
QSO must be bright enough to obtain high signal-to-noise, high-resolution spectroscopy at 
A < 4000A with a reasonable allotment of telescope time. Each one of these criteria act 
to decrease the probability that a D/H measurement can be made towards any given QSO, 
and since all the criteria must be met, the resultant probability of a QSO sightline being 
suitable for measuring D/H is very low, with only approximately 1% of QSOs at z ~ 3 able 
to provide a measurement of D/H. 

To date, there are few measurements of D/H in QSO spectra (Buries & Tytler (1998), Buries 
& Tytler (1998), O'Meara et al. (2001), Pettini & Bowen (2001), K03, Levshakov et al. 
(2002), Crighton et al. (2004)). These measurements constrain the baryon density Vt^h 2 
through the framework of standard big bang nucleosynthesis (SBBN), which predicts the 
abundances of the light elements as a function of the baryon-to-photon ratio i] and the 
expansion rate of the universe (Kolb & Turner 1990), and through the cosmic microwave 
background radiation, which provides the photon density. A measurement of the ratio of any 
of the light nuclei produced in SBBN gives the baryon density, and measurement of addi- 
tional abundance ratios test the theory (see Steigman (2006), Pettini (2006), and references 
therein for a current census of D/H and the other light element abundances). 

Recent measurements of the temperature angular power spectrum of the CMB (Spergel 
et al. 2006) also provide a measurement of Vt^h 2 depending on the assumptions made, with a 
level of accuracy roughly equal to or greater than that provided by D/H. Nevertheless, mea- 
surements of D/H are still important for a number of reasons. First, primordial D/H probes 
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the universe at one of the earliest times in the universe accessible with current observational 
and theoretical techniques. Second, the light element abundances predicted from SBBN do 
not all agree with each other; most notably the observationally inferred abundance of 7 Li is 
significantly lower than that expected from SBBN and D/H (Fields & Sarkar 2006). Third, 
D/H can help constrain deviations from SBBN, such as inhomogeneous BBN (e.g. Lara, Ka- 
jino & Mathews 2006), relic primordial particle decays (e.g. Jedamzik 2004), or non-standard 
neutrino physics (Abazajian et al. 2005, e.g.). Fourth, the dispersion in the measurements 
of D/H is larger than would be expected from the individual measurement errors (K03), i.e., 
the data demand both a better understanding of the errors on the current measurements and 
new constraints. Fifth, the value of £l b derived from D/H and SBBN requires many fewer 
priors than the CMB derived value. Moreover, the Qb from D/H can be used in principle 
as a prior in the CMB analysis, and the ratio of the values for from D/H and the CMB 
offer a precision test of the hot Big Bang model. Finally, the existing measurements of D/H 
show evidence of a trend of decreasing D/H with increasing Ahi- K03 suggested that this 
trend (and the dispersion in D/H values) is due to error under-estimation. 

Fortunately, the Sloan Digital Sky Survey (SDSS, Adelman-McCarthy et al. 2006), by 
virtue of its large sample of high redshift QSO spectra (the Data Release 4 alone contains 
5036 QSOs with z > 2.7) gives us a new data set to find those special sightlines which can 
show D/H. In this Letter, we present a new measurement of D/H in a QSO sightline from the 
SDSS, SDSS1558-0031, which was chosen as part of our high-resolution survey for Lyman 
limit absorption (O'Meara et al. 2006). 



2. OBSERVATIONS 

We have obtained two high-resolution spectra of the z = 2.83 quasar SDSS1558-0031 
using two different spectrographs, the MIKE (Bernstein et al. 2003) echelle spectrograph on 
the 6.5 meter Magellan Clay telescope at Las Campanas, and the upgraded HIRES (Vogt et 
al. 1994) spectrometer on the 10 meter Keck-I telescope on Mauna Kea. The MIKE spectrum 
was obtained as part of our ongoing high-resolution survey for Lyman Limit absorption and 
was selected from the SDSS because of the redshift and brightness of the QSO. The MIKE 
spectrum was obtained on May 10, 2004 and covers the spectral range 3221-7420 A, with an 
exposure time of 3600 seconds. The data were obtained in sub-arcsecond seeing with a one 
arcsecond slit, which provides a resolution of R = 28, 000 and R = 22, 000 on the blue and 
red arms of the spectrograph respectively. The MIKE data was reduced using the MIKE 
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reduction pipeline 1 , and has a signal to noise ratio of approximately 12 at A = 4000 A. 

The HIRES spectrum was taken on 2006 April 11 and covers the spectral range 3338- 
6200 A, with an exposure time of 4100 seconds. The data were taken in sub-arcsecond seeing 
with a 1.148 arcsecond slit, which provides a resolution of R = 34,000. The data were re- 
duced using the HIRES reduction pipeline 2 , and have a signal to noise ratio of approximately 
20 at A = 4000 A. The HIRES data are the primary source for the measurement of D/H 
presented below owing to the higher signal-to-noise and spectral resolution. With the excep- 
tion of the H I Lya line, we used the HIRES spectrum to determine all values for column 
densities presented in the text. Because we were more successful at fluxing data from the 
MIKE spectrometer, we use the flux calibrated MIKE spectrum to constrain the iVm value 
in the Lya line whose profile spans several echelle orders. 



3. Analysis 

Inspection of the MIKE spectrum of SDSS1558-0031 shows that there is a DLA at 
z = 2.70262 which is also responsible for the break in flux from the Lyman limit of the 
absorber at A ~ 3885 A. The parameters which describe the observed, single component of 
absorption (N,b,z) for the Lyman series and metal-lines in the DLA are presented in Table 1. 
The parameters and their errors are derived predominately from Voigt profile fits to the data 
using the VPFIT routine kindly provided by R. Carswell and J. Webb. For some of the 
metal-line transitions, we opt instead to use the apparent optical depth technique (Savage & 
Sembach 1991) to measure the column densities, and then use that column density as a fixed 
input parameter to VPFIT to determine the z and b values for the absorption in question. 
The absence of metal-line absorption at the position of D aruges that the observed feature 
in the Lyman series is not interploping H. 

The hydrogen column density of the DLA is sufficiently large to show damping wings in 
the Lyman a-7 transitions (Figure 1). These features allow for a precise measurement of the 
H I column density. In the case of SDSS1558-0031 the DLA Lja lies near the QSO emission 
line and the assignment of the continuum level over the full extent of the absorption feature 
is non-trivial. Fortunately, this issue is minimized in two ways. First, the MIKE spectrum is 
flux calibrated, which allows for an easier assignment of the continuum level to the Lja line, 
although it is still subject to unidentified emission features inherent to the QSO. Second, the 



1 http: / / web.mit.edu/~burles/www/MIKE/ 
2 http://www. ucolick.org/~xavier/HIRedux/index. html 
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H I derived from the damping features in the higher order lines is less susceptible to large 
continuum shape errors, because the profile spans a significantly smaller wavelength range 
than the Lyct line. In particular, the Ly/3 line of the DLA places an excellent constraint on 
the H I column density, because the line has prominent damping features, covers only m 35 
A, and has little interloping hydrogen absorption. 

To arrive at the best estimate of A H i we simultaneously vary the values of 7V H i along 
with the shape and amplitude of the local continuum level. This variation continues until 
we arrive at a value of the Nm which best reproduces the data whilst having a reasonable 
continuum shape. We adopt a redshift and velocity width of the H I, z = 2. 702646±0. 000010 
and b = 13.56 ± 1.0, from a fit to the the higher order Lyman series transitions. Of some 
concern is the fact that the redshift of the H I agrees only at the « 3 km s -1 level with the 
redshift inferred from low- ion metal-lines (Table 1). We note, however, that there exists some 
degeneracy between b the z for the Lyman series transitions we use, along with the increasing 
effects of poor signal-to-noise for shorter wavelength data (i.e. higher up the Lyman series). 
Furthermore, we cannot discount the possibility that the H I gas is multi-component, however 
there is little evidence from the metal-line transitions that this is the case. Furthermore the 
Lyman series lines all appear to be well fit using a single component, with a few departures 
due to interloping H I gas at different redshifts from the system which shows deuterium. 
When we consider the Lyman a-7 transitions, we arrive at a best estimate of the H I of 
log A ffl = 20.67 ± 0.05 cm -2 . The errors on A ffl are dominated by continuum uncertainties 
and by signal-to-noise, two effects which correlate, particularly on smaller wavelength scales. 

As can be seen in Figure 1, we observe resolved absorption by deuterium in the Lyman 
series from Ly7 all the way through to Lyman-13. Because the D I column density is large, the 
absorption is saturated until we reach deuterium Lyman-11, which offers the best constraint 
on the A DI value. For this transition, we measure logA DI = 16.19 ± 0.04 cur 2 , where 
the errors come from the error estimate of VPFIT alone (i.e. independent of continuum 
error). The transition suffers from mild contamination by interloping hydrogen on R3 25 
km s" 1 to the red side of the absorption profile, but this absorption has little effect on the 
Nj)i value. Fits to the data including and excluding a model for the interloping hydrogen 
improve the x 2 for the fit without changing the value or uncertainty in the Nm value. We 
have also estimated A DI using the AODM technique, and arrive at a consistent value log A DI 
= 16.20 ± 0.04cm~ 2 . The optical depth of this absorption feature is ideal for measuring a 
column density because it is highly insensitive to the local continuum level placement. If 
we vary the amplitude of the continuum level by as much as 20% about the adopted value 
the central value, of iV DI changes by less than the statistical error. The Nm value is further 
constrained by other Lyman series lines, e.g. the depth of the deuterium Lyman-8 transition 
rules out significantly larger or smaller values of Adi- 
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We obtain a value of z — 2.702626 ± 0.000007 for the deuterium absorption, consistent 
with that of the H I and other metals. We measure a velocity width of b = 10.48 ± 0.78 
km s" 1 for the deuterium absorption. Neutral hydrogen gas with log Ahi — 16.2 cm -2 is not 
expected to have such a narrow Doppler parameter (Kirkman & Tytler 1997) whereas the 
value is reasonable for D. 

We detect over 30 metal-line transitions in the DLA absorber. A subset of these are 
summarized in Table 1 and are shown in Figure 2. In particular, we note the presence of 
O I absorption at z — 2.702610 ± 0.000005, which is well described by a single component. 
O I absorption is important for measuring D/H in that O I directly traces the H I gas (see 
O'Meara et al. 2001), and because O I/H I O/H in most environments (especially DLA). 
Adopting the measured value of \ogN Q i = 15.86, we establish a metallicity of [0/H] = 
— 1.49 for the absorber assuming the solar (atmospheric) oxygen abundance reported by 
Grevesse, Asplund & Sauval (2005). This metallicity is higher than all the other extragalactic 
measurements of D/H. Nevertheless, a 3% solar metallicity implies minimal astration of D 
and we believe this system is still representative of primordial gas (see Figure 20 of K03 
and Romano et al. (2006)). Finally, we note that the velocity structure of the absorber, 
as traced by the metal lines, is amongst the simplest yet observed for a DLA (Prochaska & 
Wolfe 2001). 



4. Discussion 

We now discuss the value of D/H we obtain for the absorber and place it within the 
context of the combined D/H ratio for all QSO absorption systems. The best estimates 
of the H I and D I column densities in the DLA towards SDSS 1558-0031 imply a value of 
log D/H = —4.48 ± 0.06. The errors on D/H stem primarily from the effect of continuum 
placement uncertainty on the H I column density, and the signal-to-noise ratio of the data 
at Lyman- 11 where the D I column density is best constrained. 

Turning now to the combined D/H value from QSO sightlines, Figure 3 shows the new 
value of D/H from SDSS1558-0031 along with the previous values of D/H taken from the 
sample discussed in K03. We do not include the result of Crighton et al. (2004), since 
we feel that the errors on D/H in this system have been under-estimated, particularly for 
the reported N m value. We do not include the results of Levshakov et al. (2002) for the 
reasons given in K03. The horizontal solid line shows the value for the weighted mean 
of the data log D/H = -4.54581 ± 0.03606, which we round to log D/H = -4.55 ± 0.04 
to keep consistent with the literature, and the dashed lines show the ±ler uncertainties 
estimated from a jackknife analysis of the weighted means. In the case of asymmetric errors 
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on individual D/H measurements, we have adopted the larger of the errors for the calculation 
of the weighting. 

Prior to the addition of SDSS1558-0031 to the sample of D/H measurements of K03, 
a ^-minimizing linear fit to the data of the form log D/H = -2.914 ± 0.467 - (0.087 ± 
0.025) xlogiVni provided an acceptable fit to the data (P( x 2 >x | ) = 0.74). With the inclusion 
of SDSS1558-0031, however, we see a significant decrease in the likelihood that there is a 
D/H trend with AT HI . Although the data are best fit with a non-zero slope, log D/H = 
-3.707 ± 0.385 - (0.044 ± 0.021) xlogiVm, the slope differs from zero at only the 2a level 
(even before accounting for the presence of Am in both axes). Furthermore, this is not a 
good model of the data, with P( x 2 >x | ^ = 0.04. As such, the data give little confidence to 
the existence of a trend of D/H with N m . Finally, if we were to include the Crighton et 
al. (2004) and Levshakov et al. (2002) results, the likelihood that D/H depends on iVfii is 
further diminished. 

The presumption of a single value for D/H, however, is still not supported by the obser- 
vations; the observed scatter exceeds that expected assuming the error estimates reported 
in the literature. Our adopted error on the best estimate for D/H from the jackknife esti- 
mation exceeds the error on the weighted mean by a factor of two. Likewise, the x 2 of the 
six measurements of D/H about the weighted mean value of log D/H = —4.55 ±0.04 is high, 
with P(x*>Xn H ) = 0-01- Although there is the possibility that the scatter in the individual 
D/H measurements is real, we prefer the hypothesis of K03 that the errors in some of the 
individual measurements, if not all of them, are underestimated. It is likely that a combi- 
nation of new methods of error analysis and new QSO sightlines is required to fully address 
the excess scatter in the D/H measurements. 

Using the framework from SBBN (Buries, Nollett & Turner 2001), the value for D/H 
of log D/H = —4.55 ± 0.04 translates to a value for the cosmological baryon density of 
Qbh 2 = 0.0213 ± 0.0013 ± 0.0004, where the first error term comes from the errors on D/H 
explained above, and the second term from the uncertainties in the nuclear reaction rates. 
By comparison, the WMAP three year result provides an estimate of Vl b h 2 = 0.0223+q;oqo9, 
which lies within the la error estimate on VL b h 2 from D/H (Spergel et al. 2006). 

Finally, we note that the absorber showing D/H in the spectrum of SDSS1558-0031 was 
discovered serendipitously as part of our survey for Lyman limit absorption, and is the first 
D/H measurement from a QSO first discovered by the SDSS. The SDSS Data Release 3 
alone has 405 DLA with redshifts optimal for D/H (2.51 < z < 4.0 Prochaska, Herbert-Fort 
& Wolfe 2005), Assuming a small fraction of these DLA provide measurements of D/H, the 
SDSS will give many tens of measurements. The situation improves further if one considers 
the Lyman limit systems toward the SDSS QSO sample. This contrasts with the SLLS and 
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DLA which give the iVni from the Lya and Ly/5 lines, and the 7V D i from the unsaturated D I 
lines in the Lyman series. Because of this effect, the DLA and SLLS offer more pathlength 
per QSO to potentially find D/H. All of these effects combine to give a likely distribution 
of D/H measurements which is roughly independent of Nm, a hypothesis which is already 
being hinted at in the current sample, since two measurements come from LLS, two from 
SLLS, and two from DLA. Altogether, the SDSS offers the best opportunity for investigating 
the larger than expected scatter in D/H and correlations with N m , metallicity, etc. 
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Table I. IONS OBSERVED IN THE z = 2.70262 DLA TOWARDS SDSS1558-0031 



Ion 


log N a 


b 


z 


[X/H] b 


Transition A r 




cm 


km s _1 






A 



H I 


20.67 c ±0.05 


13.56 ±1.00 


2. 


.702646 


±0. 


000010 






1215.67, 1025.72, 972.54, 949.74 


D I 


16.19 c ±0.04 


10.48 ±0.78 


2, 


.702626 


±0. 


000007 






917.88 


C II 


> 14.43 d 


9.69 ±0.21 


2. 


702611 


±0. 


,000003 


> - 


2.63 


1334.53, 1036.34 


C III 


> 13.8 




2. 


.702671 


±0. 


.000054 






977.02 


C IV 


13.05 ±0.04 


9.00 ± 1.35 


2. 


702614 


±0 


.000010 






1548.20 


N I 


14.46 ±0.02 


9.20 ±0.60 


2. 


702601 


±0. 


.000005 


-2. 


.18 


1134.41, 1134.17 


N II 


13.71 ±0.04 


9.56 ± 1.13 


2. 


702660 


±0 


.000008 






1083.99 


I 


15.86 ±0.10 


6.38 ±0.40 


2. 


702610 


±0. 


.000005 


-1 


.49 


971.74, 950.88 


Al II 


12.73 ±0.05 


6.87 ±0.55 


2. 


702613 


±0. 


.000005 


-2. 


.43 


1670.79 


Si II 


14.24 d ±0.01 


7.18 ±0.14 


2, 


.702604 


±0. 


000002 


-1 


.99 


1526.71, 1304.37 


Si III 


> 13.28 




2. 


.702672 


±0. 


.000003 






1206.50 


Si IV 


12.72 ±0.04 


8.20 ±0.97 


2. 


.702595 


±0. 


,000008 






1393.76 


S II 


14.07 ±0.02 


7.05 ±0.47 


2. 


.702626 


±0. 


.000003 


-1 


.74 


1259.52 


Fe II 


14.11 ±0.03 


6.47 ±0.32 


2. 


702605 


±0. 


.000002 


-2. 


.06 


1608.45, 1144.94 



a Unless stated otherwise, the errors on the column density are from VPFIT alone. 

b [X/H] = log(X/HI) - log(X/H) where we have considered only low-ion transitions and 
have not adopted ionization or depletion corrections. 

c See text for detailed description of the values and errors for this ion. 

d Column density measurement from the apparent optical depth method. 
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Fig. 1. — H I and D I Lyman series absorption in the z = 2.70262 DLA towards SDSS1558- 
0031. The Lya transition (top panel) comes from our MIKE spectrum of SDSS1558-0031 
and the remaining panels come from our HIRES data. For the Lya and Ly/3 transitions, we 
present unnormalized data, along with a dashed line which shows our estimate for the local 
continuum level. The remaining Lyman series transitions are shown continuum normalized. 
The solid green line shows the best single-component fit to the D I and H I absorption. We 
derive our estimates of the Nm from the damping wings present in Lyman a-5, and the A^di 
from the unsaturated D I Lyman- 11 transition (~918 A rest wavelength). 
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Fig. 2. — Selected metal lines in the DLA towards SDSS1558-0031. The zero-point velocity 
corresponds to a redshift of z = 2.70262. 
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Fig. 3.— Values of the D/H ratio vs. logN m . With the new value for D/H from SDSS1558- 
0031, a linear trend of D/H with iVni is no longer statistically significant nor a good descrip- 
tion of the data. The horizontal lines represent the weighted mean and jackknife errors of 
the 6 measurements, log D/H = —4.55 ± 0.04. The right hand axis shows how the values of 
D/H translate into values for fl b h 2 using SBBN. 



